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A symbolic language allowing to solve statistical problems 
for the systems with nonabelian braid-like topology in 2+1 
dimensions is developed. The approach is based on the sim- 
ilarity between growing braid and "heap of colored pieces". 
As an application, the problem of a vortex glass transition in 
high-T c superconductors is re-examined on microscopic level. 

Statistics of ensembles of uncrossible linear objects 
with topological constraints has very broad application 
area ranging from problems of self-diffusion of directed 
polymer chains in flows and nematic-like textures to dy- 
namical and topological aspects of vortex glasses in high 
temperature superconductors In this letter we pro- 
pose a microscopic approach to a diffusive dynamics of 
entangled uncrossible lines of arbitrary physical nature. 

It is well known that the main difficulties in the statis- 
tical topology of linear uncrossible objects are due to two 
facts: (a) the topological constraints are non-local and 
(b) different entanglements do not commute. The point 
(a) can be overcomed by introducing the corresponding 
(abelian) Gauss-like topological invariant which properly 
counts windings of one chain around the other, while the 
circumstance (b) since now creates the major problem in 
constructive approach to topological theories beyond the 
abelian approximation. In order to have representative 
and physically clear image for the system of fluctuating 
lines with nonabelian topology we formulate the model 
in terms of entangled Brownian trajectories: such rep- 
resentation serves also as a geometrically clear image of 
Wilson loops in (2+l)D nonabelian field-theoretic path 
integral formalism. 

We are aimed to develop a symbolic language which 
would permit us to construct the objects with a braid- 
like topology in 2+1 dimension and to solve the simplest 
statistical problems where the noncommutative character 
of topological constraints is properly taken into account. 
The results are applied to re-examination of the problem 
of a vortex glass transition in high-T c superconductors 
P). Let us remind briefly that in CuC^-based high-T c 
superconductors in fields less than H C 2 there exists a re- 
gion where the Abrikosov flux lattice is molten, but the 
sample of the supereconductor demonstrates the absence 
of the conductivity. This effect is explained by highly 
entangled state of flux lines due to their topological con- 
straints 0. 

The most attention in our investigation is paid to a 
quantitative estimation of a characteristic time of self- 
disentanglement of a particular "test" chain in a bunch 
of braided directed chains. We distinguish between two 



situations: (i) configurations of all lines in a bunch are 
quenched and form a lattice except one test line randomly 
entangled with the others, and (ii) no chain in a bunch 
of braided lines is fixed and any chain winds randomly 
like a test one. We compute the characteristic times of 
topological relaxations r s j and r co in cases (i) and (ii) 
and demonstrate the absence of qualitative difference be- 
tween t s ; and t co , however the quantitative distinction 
between r s j and r co is shown to be sufficiently strong. 
Our result contradicts in details with the statement of 
Q on qualitative difference between r s j and r co obtained 
in the frameworks of a scaling analysis. Nevertheless our 
result does not destroy the physical conclusions of the 
paper || about the possibility of topological glass transi- 
tion in entangled flux state in high-T c superconductors. 
According to the above mentioned cases (i) and (ii) we 
define respectively two discrete models I and II. 

The model I is as follows. Take a square lattice in 
(sy)-plane with a spacing c and put in all vertices of this 
lattice the uncrossible obstacles. Consider a symmetric 
random walk with the step length c on a dual lattice 
shifted by ~ in both x- and y-directions. We are in- 
terested in computing the probability P S i{N) of the fact 
that after N steps on the dual lattice the random path 
will be closed and unentangled with respect to the ob- 
stacles. It is clear that in the (2+l)-dimensional "space- 
time" 1. 2 (x,y) x Z + (i) this model describes statistics of 
a "world lines" (time-ordered paths) of a single particle 
jumping on a square lattice in (x7/)-projection, making 
each time a step toward i-axis and topologically interact- 
ing with the lattice of infinitely long straight lines. The 
(xy)-section of this model is shown in fig.fy. 



FIG. 1. The (a;j/)-projection of a path unentangled with 
the square lattice of topological obstacles. 

It is well known that there exists a bijection between a 
path in the lattice of obstacles with coordinational num- 
ber z and a path on a Cayley tree with z branches . Dif- 
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ferent topological states of a given path coincide with the 
elements of the homotopy group of the multi-punctured 
plane T^, generated by a countable set of elements. The 
translational invariance allows to consider in a local basis 
the group Tqo/Z 2 = T z / 2 , where T z / 2 is a free group with 
z/2 generators i.e. a z-branching tree. Any topological 
state of a path in the lattice of obstacles is encoded by an 
element of the group T z / 2 , and uniquely corresponds to 
some irreducible word written in terms of generators of 
the group T z / 2 . The length of this irreducible word coin- 
cides with the geodesic distance along the Cayley graph 
of the group T z / 2 . The average "degree of entanglement" 
of an A-step path in the lattice of obstacles with co- 
ordinational number z for A 3> 1 is characterized by 
the averaged geodesic distance || (L s i(N)) = ^jr-N on 
the Cayley tree with z branches. Hence, the normalized 
"complexity" (Z S i) of a typical topological state of a path 
can be defined as: 
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For diffusion on the Cayley tree, P s i(fc, N) deftermines 
the probability for the A-step symmetric random walk to 
have a distance between ends in k steps along the tree. 
This probability has been computed many times — see, 
for example, jl. Thus we reproduce the final result for 
the return probability P s i(k — 0, A) as A — > oo on the 
Cayley tree: 



P si (A) = P ai (k = 0,N) 
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where q = - and p = 1 — q = ^— ^ are the probabilities 
of steps "towards" and "backwards" the origin of the 
Cayley tree; a = 2^/pq. Define A, a span of the A-step 
random walk in (a;?/)-projection. In physical terms of 
the original " vortex problem" , VA 2 " is the averaged size of 
thermal fluctuations of the vortex line. Hence, we can set 
N = At and estimate the time of topological relaxation 
T s i = n run of a single vortex line in an ensemble of 
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In contrast to the model I, the model //describes collec- 
tive dynamics of the world lines and ultimately leads to 
the consideration of the (2+l)~dimensional ("surface") 
braid group P 2 +i- Consider the two-dimensional lat- 
(x, y) and take distinct points Pi, P 2 , ...P( n+ i)2 , e 
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A (2+l)-braid of (n + 1) strings on Z based at 



{Pi, .... P( n+ i)2} is an (n+ l) 2 -tuple b 
of paths, such that: (1) &t(l) = 
{Pi,...,P„ 2 } Vz G + l) 2 } 

V{i,i}e{l, 



Pi and 
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(2) h{t) ± VJ 
(n + l) 2 }, (i f^j, te [1, A]). The braid 
Z 2 based at {Pi, P(„ +1 )2} is the group 

-P(n+1) 2 }- 



group P„+! on 

of homotopy classes of braids based at {Pi 
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The group 

{l,...,n+ l},j 
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has 2n 2 + 2n generators a. 
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/a + 1} and their inverses with the standard rela- 
The geometric representation of generatoes of 
is shown in fig.pja. An element of the braid group 
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is set by a word in the alphabet {o~\ x ,a 
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By the length A of a braid we call a length of a word in 
a given record of the braid, and by the irreducible length 
(or primitive length) — the minimal length of a word, in 
which the braid can be written. The irreducible length 
can be also viewed as a distance from the unity on the 
graph of the group. Graphically the braid is represented 
by a set of strings, going upwards in accordance with a 
growth of a braid length — see fig.||. 
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FIG. 2. The generators of the surface groups B^ D (a) and 
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FIG. 3. The (2+l)-dimensional braid. 

Define now a symmetric random walk on a set of gen- 
erators ("letters") {cr^, <J\i , ■ ■ ■} with the transitional 
probability 2 n' i +2n • Namely, we raise recursively an 
A-letter random word W adding step-by-step the let- 
ters (say, from the right hand side) to a growing word. 
The probability that the A-letter word is completely 
contractible (i.e. has a zero's primitive length) defines 
the probability to have topologically trivial braid of the 
record length A. 

Our main tool in the investigation of the braid group 
is the so-called locally free group ||. The (2+l)D ("sur- 
face" ) locally free group CT 2 ^ is obtained from the braid 
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group Bl» x by omitting the braiding relations. The * (h co (N\CT 2 n D )) < (L^N^ 2 ^)) < (L co (N\CT 2 n D )) 



group CT 2 ^ has 2n 2 + 2n generators /y , i G { 1 

l},je {l,...,n} ; /# ) ,ie{l,...,n},je{l,... 
and their inverses with the relations (see also fig.l 
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There is a bijection between words in locally free group 
and colored heaps, whose elements are either "white" 

. That is, any word written 
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/£ ■»> or "black" (./;;•* 

in terms of letters-generators of the group CT^ repre- 
sents a configuration of a colored heap (see fig.|J) in a 
box of the base ofnxn cells and any such heap uniquely 
defines some word in the group CJ 7 ^ . The configuration 
of the heap with a "black" block following immediately 
after a "white" one in the same column is forbidden. 
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FIG. 4. The (2+l)-dimensional colored heap. 

For the uniform Markov dynamics on the set of gen- 
erators we can compute the average primitive length 
(L co (N)) of the iV-letter word which characterizes the 
degree of entanglement of threads in a braid (compare to 
the definition of (L s i(N))). Let us take into account that 
CT 2 ® is a subgroup of -B^+i, while B^_ x is a factor- 
group of CT 2 ^ . This is manifested in the two following 
facts. (1) By definition the locally free group CT 2 ° has 
less relations than the braid group B^+i- Hence, the 
number of distinct words of primitive length L co in braid 
group is bounded from above by the number of distinct 
words of the same primitive length L co in the locally 
free group. (2) By construction (compare figs.^a and ||b) 

f^: v) = [a] ' -' 1 j e [l,n]). Thus, the number 

of distinct words of length 2L co in the braid group is 
bounded from below by the number of distinct words of 
the length L co in the locally free group j8|. The facts 
(l)-(2) allow us to get the bilateral estimation for the 
average primitive length (L co (N \B„+i)) of the iV-step 
random walk on the surface braid group B^P_ X . 



(4) 

The computation of ( x L co {N\CJ- 2 ^')) involves the con- 
cept of the roof of the heap. In physical terms a roof of 
a heap consists of a set of " most top" blocks which can 
be removed from the heap without changing the rest of 
it. The projection of a roof to the (a;y)-plane for some 
particular configuration of the most top blocks is shown 
in fig.|[ Let us stress that local hights (measured from 
the bottom of the box) of different roof's blocks might 
be different. 



FIG. 5. A particular configuration of a roof of (2+l)D heap 
is shown by bold segments. 

The process of growth of a heap (i.e. the random walk 
on the surface locally free group CT 2 ^) consists in adding 
step- by-step new "black" or "white" blocks to the roof. 
Hence the dynamics of a heap is controlled by the dy- 
namics of a roof. For a particular configuration of a roof 
we define the "size" of a roof #T (i.e. the number of bold 
segments in fig.|^) and the number of empty segments m 
having i bold neighbors (apparently, rii = OVi > 3). In 
figJI one has n = 5, n = 6, n.1 = 25, n 2 = 17; #T = 12. 
For the values #T, rii the following conditions hold: 



n + m + n 2 
6#T- 8(n + 



1) < ni 



2n 2 + 2n 
+ 2n 2 < 6#T 
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For a given configuration the local dynamics of a size of 
a roof reads 

A#T = 1 with probability 
A#T = with probability 
A#T = -1 with probability 

which together with (|^) allows to estimate #T as follows 



1 - <A#T) < 



2n 2 + 2n 



< 1 



8(n + l) 
2n 2 + 2n 



- <A#T) 



In a stationary case we have (A#T) = 0, what permits 
to get the asymptotic value (#T) of the average number 
of block in the roof for n ^> 1 : 



(#T) 
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Return to the random walk on CJ- 2 ® ■ The conditional 
change of the primitive length L co (N\CJ :2 ®) for one step 
of the random walk is 



&L co (N\CT 2 n D ) = 1 with probability 1 



&L co (N\CT 2 n D ) = -1 with probability 



# 



2(2n 1 +2n) 



2(2n 2 +2n) 



Taking into account (^) we get in a stationary state (when 
N ^> 1 and n 3> 1) the following asymptotic value for 
the average primitive length (L co {N\CT 2 f')y. 



(L co {N\CT 2 n D )) = -N 



Thus, according to (Q) we arrive at the bilateral estima- 
tion for the average length of the primitive word for the 
TV-step random walk on the surface braid group 
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The quantity (Z c 
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characterizes 



the "complexity" of entangled state. Comparing (0) to 
(Q) we can conclude that the averaged topological state of 
a braid (l/ co (TV|_B^ D )) obtained in course of a collective 
motion of lines has the same asymptotics in N as the one 
of a single line motion and interpolates between entangled 
states in effective lattices of obstacles with coordinational 
numbers z e s ranging in the interval [|, 14]. 

Let us utilize now the "heap concept" and eq.@ for 
estimating from above the time of topological relaxation 
in a bunch of vortex lines, considered as a braid of di- 
rected random walks. The random growth of a braid 
of 2n 2 + 2n lines of length N each, can be interpreted 
as an JV— step random walk on a surface braid group 
-B^+i ■ The topological state of a braid is uniquely char- 
acterized by a primitive word Wn in terms of genera- 
tors of -B^+i- The disentangled state of two neighbor- 
ing trajectories means (in terms of the group -B^+i) that 
the primitive word Wjv does not contain the generators 

°th°ih » HJ) 1 for some 1 - {hj} - n - 

Let L,j(Wn) be the number of generators a = 

'l' 7 *' "j ' i^ij'*) j" m ^ nc P r i m iti v e word Wjy. Then the 

function P co {x, N} = P{L ? (Wn) = x} defines the char- 
acteristic time T ro = 



Pco(0,JV) 

particular line in course of collective Brownian motion 
of all 2n 2 + 2n lines in a braid. Suppose that after 
N random steps we arrive at the word Wn- Then at 
the step N + 1 we have Wn+i = Wn Wnn+i- The 
probability that the word Wnn+i contains a genera- 
tor from the set a (with a prescribed sign) is P{a £ 
Wnn+i} = \ x \ = §5- Consider the probability 
q b N = P{L a (W N+1 ) = L a {W N +i) ~ 1} of reducing the 
element a by 1 at the step N + 1. The probability q b N we 
can estimate from below, replacing the word Wn+i by 



the word Wn+i which itself is obtained by replacing the 
generators of the group B 2I ^_ 1 by the ones of the group 



CT 2 ^ . The word Wn+i is also primitive by construction 
and q l N = P{L f {W N +i) = L f {W N +i - 1)} < q b N , where 
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Denote by Q l N the probability 



that after N steps the roof contains the generators /, 



or /,, 
ary i 

3 G {1 



. In the case n » 1 we can neglect the 'bound- 
ary effects' caused by marginal generators f^ + \ • , f[ X j , 



probability Q l N is the same for all / G ^ /. 



, n} so that the 
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Thus we have the obvious expression for the average num- 
ber of blocks in the roof in the stationary case N 3> 1 : 
(#T) = 2n 2 Q l ] which gives Q l = \. Then the probabil- 
ity of cancellation of the generator / at the step N + 1 is 
In = T&Qn = 955 — Q b where q b is the stationary prob- 
ability of cancellation of generator a at N — > oo. Substi- 
tuting (7 fc for g in (||) we get the estimate from above for 
the time r"P; 
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The numerical analysis of statistics of braid and locally 
free groups || enables us to conjecture that the estimate 
r" p for the group CT 2 ^ is close to the value r co for the 
group B 2D . 

Comparing (||) to (||) we conclude that the scaling de- 
pendences of the characteristic times r S i and t co on — are 
the same, however numerically r co is much larger than r 8 i, 

2 

namely lim -p- In ^ ss 6.5. 

■7— 31 

Summarizing the said above, let us emphasize that the 
random braiding can be analyzed within the frameworks 
of a symbolic dynamics on the locally free group describ- 
ing the growth of a heap of colored pieces. The proba- 
bility to have disentangled state of two vortex lines can 
be estimated from above by the probability to have no 
pieces (blocks) in a given column of a heap. This model 
seems to be a natural discretization of a standard ballistic 
growth process of Kardar-Parisi-Zhang type || . 

The authors are grateful to A.Vershik and J.Desbois 
for useful discussions. 
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